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'5 Abstract 

Superfluids support many different types of sound waves. We investigate the relation 
between the sound waves in a relativistic and a non-relativistic superfluid by using 
hydrodynamics to calculate the various sound speeds. Then, using a particular holo- 
graphic scalar gravity realization of a strongly interacting superfluid, we compute first, 
^ second and fourth sound speeds as a function of the temperature. The relativistic low 

temperature results for second sound differ from Landau's well known prediction for 
QQ the non-relativistic, incompressible case. 
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1 Introduction 

In this work, we are interested in the sound speeds of a strongly interacting relativistic 
superfluid. Since perturbative techniques are not applicable, we use the AdS/CFT corre- 
spondence [U 12113] to study the sound modes in these systems. Relativistic superfluids may 
be important for understanding the physics of neutron stars [1]. In addition, the holographic 
model we use has also been employed to model quantum criticality in certain condensed 
matter systems (see O [6] for reviews). 

A superfluid can be defined as a liquid in which a Bose condensate is formed. In 
a weakly coupled system, a Bose condensed phase occurs whenever a large number of 
particles occupies the ground state of the free theory. A more formal description of a Bose 
condensate is that of a system where the single particle density function does not vanish at 
large spatial separation [7]. The virtue of the latter definition is that it does not rely on 
a weak-coupling description of the system. Indeed many properties of the superfluid phase 
of ^He can be accounted for when modeled by a strongly-interacting Bose condensate. Of 
particular interest to us in this work are the sound modes which may be excited in the 
superfluid. 

A superfluid at nonzero temperature can be thought of as a two component system: 
an uncondensed, normal component and a condensed, superfluid component. Each compo- 
nent has its own density field and velocity field. The collective motion of the fluid where 
both components move in phase is called ordinary sound, while second sound is associated 
with out of phase motion of the two components. In a non-relativistic and incompressible 
superfluid one finds that ordinary (or first) sound waves couple more strongly to pressure 
oscillations while second sound waves are sourced by temperature oscillations [8]. There is 
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also a fourth sound mode where the normal component is immobilized, for example by pack- 
ing a capillary tube with powder. This sound mode propagates through density fluctuations 
of the superfluid component]^ 

Regardless of the coupling strength, a hydrodynamic description of the superfluid en- 
ables a calculation of sound speeds from the equation of state and various thermodynamic 
quantities. In section [2] we describe how such a calculation can be carried out for the rel- 
ativistic superfluid and relate the relativistic expressions for the sound modes to the much 
better known non-relativistic ones. 

The particular holographic model we use is the scalar gravity system introduced by refs. 
|10^ [TT] which consists of an Abelian gauge field and a charged scalar field in an electrically 
charged black hole background with a negative cosmological constant. A Higgs mechanism 
in the bulk gravity theory is dual to spontaneous breaking of a global U(l) symmetry in 
the boundary field theory. The authors of ref. [12] argued that this system is dual to a 
strongly interacting, relativistic superfluid. We give more details regarding the mapping of 
the gravitational system to the conformal superfluid in section |3j 

Our work was motivated in part by Landau's well known prediction for the low temper- 
ature behavior of second sound in a non-relativistic incompressible superfluid |8j. Landau 
argued that at low temperatures, 

lim = ^ , (1.1) 

where Cj is the nonrelativistic speed of i'th sound and d is the number of spatial dimensions. 
This argument has been generalized to a relativistic fluid in [13^. In our setup, we find 
that the speed of second sound does not approach this limit. In section |4] we describe 
the behavior of second (and fourth) sound for the holographic superfluid. We discuss the 
difference between ( 1 . 1 1 and the relativistic result in section |5} 

Previous work [12l [HI [151 Ull E] examined the speed of second sound and fourth sound 
in a probe limit where the scalar and gauge field were not allowed to back react on the 
metric. (The probe approximation can also be thought of as a limit where the charge of 
the scalar field becomes large.) This limit is not well suited to an investigation of Landau's 
prediction. The reason is that in the probe limit the charged matter is a small perturbation 
of the rest of the system, while Landau's prediction requires an assumption that the whole 
system become a superfluid at T = 0. In the probe limit fourth sound and second sound 
coincide and approach first sound at zero temperature provided that the zero temperature 
limit of the theory is well defined. In this paper, we do not work in the probe limit and 
consider the full backreacted geometry. 



2 Sound modes 

Consider a thermal system with a conserved U{1) symmetry. The hydrodynamic variables 
used to describe such a state are the local velocity field of the fluid u'^ with = 0, . . . , d, 
the temperature T and the chemical potential /x. If the U{1) symmetry is spontaneously 
broken a condensate forms, and the resulting Goldstone boson ip provides for a new degree 

^ Third sound is associated with surface waves in a thin film of superfluid [9]. 
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of freedom = d'^ip. Refs. |18| [T2] developed a relativistic hydrodynamic description of 
such a system in the ideal limit, ignoring dissipation]^ In what follows we summarize their 
results. The extra degree of freedom is interpreted as the (non-normalized) velocity of 
the condensate v'^ via ^'^ = fiv'^. The norm of the condensate velocity, (3, defined through 
(1 — = —v^v'^ serves as the new hydrodynamic degree of freedom. The velocity field 
n*^ denotes the velocity of the uncondensed phase. It can be defined as the local boost 
parameter required to bring the superfluid velocity to the form = (l,/3n) where n is a 
spatial 3-vector of unit norm. The pressure P, being a Lorentz scalar, can depend on T, 
fi and (3. We define the entropy density s, normal density pn and condensate density ps as 
the variables conjugate to temperature, chemical potential and the norm of the superfluid 
four-velocity, 

dP = sdT + pndp - ^df , (2.2) 

where = CuC^- When gradients of the hydrodynamic variables can be neglected, the 
resulting stress tensor and U{1) current of such a system take the form 

T'^- = {e + P) u^u'' + Pr,"'' + mPs^"^^" , , , 

(2.3 

J = PnU + PsV 

where e is defined as the Legendre transform of the pressure with respect to the temperature 
and chemical potential, e = —P + Ts + p^p,. With this definition, e differs from the time- 
time component of the energy-momentum tensor. In what follows we will use w = e + P 
and p = Ps + Pn- In the superfluid literature, there exists an alternate definition of the 
condensate density via a current- current correlator. That these two definition coincide was 
shown in ref. [22]. For convenience, we reproduce the result in appendix [A| 
The hydrodynamic equations are given by the conservation equations, 

d^T"'' = , d^r = . (2.4) 

These conservation conditions are supplemented by a Josephson equation: 

u^vu = -1 . (2.5) 

One way of deriving the Josephson condition along with the conservation laws (2.4) and 



eqs. (2.3) is the Poisson bracket technique used in refs. |18|I22]. Note that with the absence 



of dissipative terms in eqs. (2.3), eqs. (2.4) and eq. (2.5 ) imply a conserved entropy current: 



d^, {u^s) = . (2.6) 



2.1 Second sound 



To look for sound modes which propagate in the system described by eqs. (2.4) and (2.5) 



one should study linear perturbations of the hydrodynamic variables around a static con- 



figuration. The system (2.4), (2.5) supports two sound modes which are given, in the ideal 



limit where damping effects are ignored, by the two positive roots of 

ac^ - /Jc^ + 7 = (2.7) 



^See refs. [T^[^[1T1[^ for earlier work on the relativistic hydrodynamics of superfluids. 
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(We follow ref. in our treatment of the sound modes.) If the system is conformal so 



that the only dimensionless quantity is p/T, then one of the two positive roots of eq. (2.7 1 
is given by 

(2.9) 



1 

d' 



The mode associated with this phase velocity is called normal sound. The other of the two 



roots of eq. (2.7) is then given by 

sTps 



4 



dp 
dp 



1 fds\ 
s \dp)^ 



-1 



(2.10) 



Wpp yp yiyp,/ y 

and is called second sound. If instead of the entropy per unit volume s we use the entropy 



per particle a = s/p^ then eq. (2.10) reduces to 



1 



w 



{da/dT), 



(2.11) 



where we have used {da/dp)T = —{T/p){d(j/dT)^ which follows from scale invariance. 

The relativistic formulation of hydrodynamics developed in ref. [18] and described above 
coincides with the relativistic version of the Landau-Tisza two fluid model for superfluids 
\2U\ 1^ . In what follows we would like to compare eq. (2.11) with the non-relativistic 



expression for second sound obtained directly from the Landau-Tisza two fluid model 
(see for example ref. |26] for a review). In the non relativistic setup first and second sound 
are given by the solutions to eq. (|2.7|) with 



a 
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1 



1 



m 



1 Ps 



dP^ 
dp 

2 fdT\ fdP\ 



dT 
d^ 



(2.12) 



m'^ Pn'^ [^da ) Adp ) . 



where m is the mass of the condensing particles and p, ps and pn are the total number 
density, superfluid number density, and normal component number density respectively. 
The quadratic equation is typically solved by making an additional assumption about the 
compressibility of the fluid, kt = where Ka = {dP/dp)a and kt = {dP/dp)T The 
difference kt—Ku is proportional to the square of the thermal compressibility of the material. 
For superfluid ^He, the thermal compressibility is tiny (1 — k^/kt ~ 10~^), justifying this 
assumption and the use of the term "incompressible". With this assumption, the non- 
relativistic result for second sound is 



1 



mp^ {da/dT), 



(2.13) 
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We wish to compare the incompressible result, eq. (2.13), with the scale invariant result. 



eq. (2.11 ). To that end, we define the non-relativistic limit by taking the chemical potential 
to be very large and approximately equal to the mass of the particles, and the number 
densities approximately equal to the corresponding charge densities: 



m 



w 



Pn ~ Pn 



Ps ~ Ps 



(2.14) 



In this limit, indeed the scale invariant result (2.11) is close to the incompressible, non 



relativistic result. The analysis of ref. p7] anticipated such a connection. 



By directly comparing the coefficients (2.12 ) and (2.8 ), we can recover the non-relativistic 



sound speeds from the relativistic result in a more general context that does not assume 
scale invariance or incompressibility. By replacing the dependent variables p and T in the 



coefficients (2.8) with a and p, we find, not assuming scale invariance, that 

fdT\ 



a = 1 , 

Pn 



f3 



w 



7 



dP 



p w \d(T J 



dp/T 



Ps crT\ 
» I \^ + ——) -2 

9pJa\ Pn P J 



dP\ p^a^T 



+ 



dT 



2Ps 



a 



aT 
P 



(2.15) 



Comparing the coefficients (2.12) with the coefficients (2.15), the corresponding quadratic 



equations (2.7) agree in the norelativistic limit (2.14). 

For the numerical work in section |4j instead of working with a and p as hydrodynamic 
variables, we find it convenient to use the chemical potential and entropy. To convert from 



the T and p variables of (2.8) to the s and p variables, we note that in the scale invariant 
case 



dp^ 


\ pd - pp' 


(dp 


df^ 


K T-T'p' 


\dp 


ds^ 


\ sd 


fds 




l^-T-T'p ' 


[dp 



Tp' - T'pd 
T-T'p 

sT'd 
~ T-T'p ' 



(2.16) 



where T' = {dT/dp)s and p' = {dp/dp)s- Equality of mixed partial derivatives implies that 
{ds/dp)T = (dp/dT)^ and hence that 



— sT'd + pp = pd . 



Given these expressions, it is easy to see that 

SPs I 



wpp 



T - T'p) 



(2.17) 



(2.18) 



Near the critical temperature, second sound vanishes due to the vanishing of the superfluid 
density ps- 

In non-relativistic superfluid ^He, it is expected that at low temperatures second sound 
will approach the speed of sound in a phonon gasj^ This expectation is based on the 



^If the sample is not pure, i.e. ps does not approach p as T ^ 0, the second sound speed should vanish 
at low temperature [211 ESI ■ 
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assumption that at low temperatures the dominant excitations are phonons whose dispersion 
relation is given by u; = kcq. In ref. [13] it was argued that in a relativistic setting such a 
dispersion relation will lead to 



sT 



sT + ;U/Jn 



(2.19) 



and to T {ds/dT) = sd. Implicit in the identification (2.19) is the assumption that only 



the phonons contribute to s and pn at low temperature. Using eqs. (2.16) and (2.18), and 
assuming that p ~ /x"^ at low temperature, we find 



lim Co 

T~>o ^ 



Z3_ 
d 



(2.20) 



This expression remains unchanged in a non-relativistic setting. See for example ref. |30] . 
2.2 Fourth sound 



If the normal component of the fiuid is prevented from moving, it is possible to excite a 
different kind of sound mode, fourth sound [31] |9]. To observe fourth sound experimentally, 
superfiuid is channeled through a tube packed with a powder that immobilizes the normal 
component [32]. In such a setup momentum is not conserved; when computing the sound 
velocity we need, once again, to find the phase velocity for linearized fiuctuations around 



a static background, but here omit the momentum conservation equation (2.4). In the 



formulation we are using, ref. |17j carried out such an analysis and found that the phase 
velocity associated with fourth sound is 



2.21 



Ref. [12j observed that fourth sound asymptotes to second sound in a large temperature 



limit. To make this relation more transparent, we use eq. (2.21 ) to rewrite eq. (2.11 ) in the 
form 



cl 



C2- 



W 



ds_\ 
df), 



At high temperatures, we expect s ^ T"^ and w ~ sT. 

The non-relativistic expression for fourth sound, computed in ref. [9 , is given by 



'2.22) 



Ps ~2 I ~2 
— Cl + —02- 

P P 



(2.23) 



To compare eq. (2.23) with eq. (2.21 ), we use eq. (2.17) and eq. (2.18) to obtain p' in terms 
of c^. We can then rewrite eq. (2.21) as 

PPs 



cl 



-cl + 



w 



sT + pp sT + pp 



(2.24) 



In the non relativistic limit (2.14), the expression (2.24) reduces to eq. (2.23). 

Ref. [32] observed experimentally that in ^He, fourth sound closely follows eq. (2.23): 



C4 vanishes at the critical temperature, and approaches first sound at low temperatures. A 



similar behavior is expected from the relativistic formula (2.24). At low temperature and 
for a pure sample for which ~ p we obtain C4 ~ cf . 
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3 The Gravity Dual of a superfluid 



In this section, with an aim of calculating sound speeds in mind, we consider a holographic 
model of a strongly interacting, relativistic, scale invariant superfluid [lOl [11] [12]. The 
action for a Maxwell field and a charged complex scalar field coupled to gravity is 



S — S'bulk + 'S'boundary (3.25) 



where 

-Sbuik = I d^x 
with 



(3.26) 



mV'l) = m2|V^|2 + ^|V^|4 (3.27) 

and ^boundary IS a boundary action which ensures a well posed variational problem and also 
renders the on-shell action finite. We reveal the various terms in S'boundary below as we 
need them. Unlike a traditional four dimensional Landau-Ginzburg model, a u 7^ is not 
necessary to see a superfluid phase transition in this gravitational system. However, since 
the zero temperature description of the model with u ^ seems to be under good control 
|33j , we will be interested in seeing the effect of a higher order term on the low temperature 
behavior of the systemj^ Roman indices a,b, . . . run from to 4 and are raised and lowered 
with the five dimensional metric gab- The Greek indices ^, z^, . . . used in the discussion of 
sound speeds, and which run from to 3, are raised and lowered with the Minkowski tensor 
r]^'^ with signature ( — h ++). 

We look for solutions to the equations of motion that follow from this action that are 
asymptotically anti-de Sitter with a flat slicing. More specifically, we require the metric to 
approach the form 

ds^ = ^{-df + dxj) + (3.28) 

at large r. The conformal boundary of the space is the constant r slice at r — > 00. From 
the usual rules of the AdS / CFT duality [2] 13] the boundary value of the metric g^y acts as 
a source for the boundary theory stress tensor T^^; the bulk gauge field acts as a source 
for a conserved current corresponding to a global U(l) symmetry; and the near boundary 
data of the scalar ip sources a scalar operator Oa with conformal scaling dimension A > 1 . 
The conformal dimension of Oa is related to the five dimensional mass of the scalar field 
m? through m?'Lp' = A(A — 4). In what follows, we will be considering the particular case 
of m? = —15/4^2. This mass yields operators of dimension of either 3/2 or 5/2 depending 
on which term in the near boundary series expansion of ijj we choose to be our source |35] . 
By boundary values, we mean more specifically 

gfl = hm {L/rfg^, , Af = lim A^ , = lim {r / L^f'^^i^l^L^'^) . (3.29) 



''in [311 a different bulk action was considered where a global (7(1) symmetry is broken via the Stuckelberg 
mechanism. It would be interesting to study the behavior of second sound in such a setup. 
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Given these definitions, the one-point functions of the corresponding field theory operators 
are 

;^/,,/v , 2 5S , , 1 5S ,^ V , 1 SS 

™ = hm —TV , = lim , y-r , (Oa) = lim , y-r . 

^ ' r^^y^6gl^y r^^y^S^ib) 

(3.30) 

We will see shortly that {T^'^), (J^'), and (Oa) are proportional to the coefficients of the 
0(r^^), 0{r~'^), and 0{r^'^) terms in a near boundary series expansion of the metric g^y, 
gauge field A^, and scalar field respectively. 

The utility of this action is that it describes a system that undergoes a superfluid phase 
transition. If the bulk scalar field vanishes then the equations of motion following from eq. 



(3.26) admit a Reissner-Nordstrom black hole solution. Such a gravity solution corresponds 
to a thermal phase of the field theory with a nonzero chemical potential. Following the 
general arguments in ref. PIJ] once the temperature of the black hole falls below a certain 
critical value Tc, a new phase exists in which the scalar field condenses. A solution with 
a nontrivial profile for ip in the bulk corresponds to a superfluid phase in the field theory. 
In what follows we will construct such solutions and map bulk to boundary quantities 
explicitly. We first describe the construction of a static and stationary configuration of the 
superfluid which is similar to the solutions constructed in refs. \36\ 157] . In section 3.2 we 
study linear vector perturbations of this configuration which will enable us to obtain the 
remaining thermodynamic quantities required to compute the various sound modes. 

3.1 Static and isotropic configurations 

First, we construct static and isotropic configurations. Such configurations are described 
by the metric and gauge field 

ds^ = -f(r)dt^ + — - + — rdxf , A = <p(r)dt (3.31) 
g{r) 

where /(r) and g{r) have a simple zero at r = tq which is a black hole horizon. We make 
the gauge choice that ip{r) is a real function. The equations of motion for ip, (p, /, and g 
are 

^ ^ \2g^ 2f^ rj ^ ^ fg ^ 2g ' 

u, I a' f ^\ u 2eV^^ 
= (/)" + - ^ + * ^' ^ ^ 



2g 2f rj^ g ^ ' 

.n'^^^Jl 3/^^gW 

^44rg 44rf^ fg ' 

((p'f 3g' 3f' V 6 3 



(3.32) 



Near the boundary we require that / = r"^ / + 0{r^), (f) = ii + 0{r^^) with /x a constant 
and that the coefficient of the 0{r^~'^) term in a series expansion of ip vanishes. The latter 
condition ensures that the boundary theory scalar field is not sourced. These boundary 
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conditions imply the near boundary series expansions 



r\ 3/2 / j2 
-) + + 

r J \ r 



(j) = 



"2^ 



fir) 
9{r) 



r 

1? 
^2 



(3.33) 



+ 



r2 KiL(V'(°))2 2/^2l3P-5k2l3^(0)^(i)/2 



+ 



+ ... 



L'^ r 

with either il)^'^'^ or '(/^(^^ set to zero. The coefficients P, and V'^"'* or ^/^^-^^ are integration 
constants which are determined by requiring that /, g and (p vanish at the black hole 
horizon located at r = rg, and that '(/' is finite there. Since the boundary value of the gauge 
field sources the U{1) current, we may interpret the boundary value of At as the chemical 
potential ^. We have used the letters P and p in anticipation that these are the pressure 
and charge density respectively. 

The partition function Z for this configuration can be computed by evaluating the on- 
shell Euclidean action The Euclidean action is different from the Lorentzian action 



given in eq. (3.25) by a minus sign: 



InZ 



-Se = S 



(3.34) 



The bulk on-shell action (3.261 reduces to a total derivative. Introducing a UV regulator 
Too we find, 



'S'bulk 



1 



2k2L3 



(3.35) 



The boundary action in (3.25) is composed of a Gibbons-Hawking term which ensures a 



well posed variational problem for the metric 

1 



^GH 



24 



d'^x^J^{2K -Q/L) 



(3.36) 



In eq. (3.36), 7 is the induced metric on the large radius slice r = Too and K = j Vani, 
is the trace of the extrinsic curvature, (n" is the outward pointing unit normal vector to 
the boundary.) Additional terms in ^boundary are required in order for the boundary terms 
associated with the variation of the scalar field to vanish on shell. Thus, we have 



5. 



3/2 



'5/2 



d^x 



-7 ( 2V'n'^9aV' + li^VL 



-7 



3 ,2 



(3.37) 



Using the boundary action 



•^boundary — SqH ~\~ 



(3.38) 



^ Strictly speaking, Euclidean time should be compactified on a circle with a radius equal to the inverse 
temperature. For brevity, we avoid writing this out explicitly. 
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we find that the partition function is 



\nZ = P 



Volume 



(3.39) 



where we have assumed that either i/^^^) or ^/^^^^ have been set to zero. As anticipated, the 



integration constant P in (3.33) may be interpreted as the pressure. With the boundary 



action (3.38 ) we can use standard techniques [21 El EHl [33 SOI SI] to compute the one-point 



functions (3.30). We find that 



(T^'^) = diagonal {?,P P P P^ , ( J^) = (^p o) , 

(O3/2) = -2V^°^ , (O5/2) = 2^^'^ . 



(3.40) 



The details of this computation are left to the appendix. Apart from the pressure, chemical 
potential and charge density, other useful thermodynamic quantities are the temperature 
and entropy. The entropy of the configuration is given by the entropy of the black hole 
which reads 

° (3.41) 



In a similar manner, the temperature is given by the Hawking temperature: 



47r 



r=ro 



(3.42) 



The absolute value of the expectation value of the scalar field |(Oa)| can be interpreted 
as the order parameter for the phase transition [TT]. Thus, with a static and stationary 
configuration, we are able to determine all the thermodynamic parameters of the superfluid 



except for the superfluid density. We obtained explicit solutions to eqs. (3.32) numerically 
and discuss some details of the results in section HI 

3.2 Vector perturbations of the static and isotropic system 

In order to determine the superfluid density ps, we consider a configuration where the 
superfluid velocity Vi, and normal component velocity Ui, are small, unequal and in the 'x' 
direction (we define x = (x, y, z)). In such a setup the space-time component of the energy 
momentum tensor and the space component of the gauge field are related to Ux , Vx and ps 
through 



(e + P)Ux + p.PsVx = {sT + ppn)ux + pPsVx , 



= p^Ux + PsVx 



(3.43a) 
(3.43b) 



As described in section [2] Vx is proportional to the gradient of the Goldstone boson 93: 

Cx = fJ-Vx = dx^p- (3.44) 



Thus, by computing {T ), (J^) and S,x, we can use eqs. (3.43a), (3.43b) and (3.44) to obtain 



Vx, Ux and, in particular, pg. To generate a solution where Ux and Vx are small, we consider 
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linear fluctuations of Ax and Gtx around the background (3.31 ). Instead of the background 



metric and gauge field in (3.31) we use 

^2 



-f{r)dt' + 



2r' 



I ^ + T9'^^i + -prGtxdtdx 
g[r) 



A 



{r)dt + Ax{r)dx (3.45) 



and work to linear order in Gfr and A. 



The linearized equations of motion for Gtx and Ax around the background (3.45) are: 

2eW, ^r^G[x 

n f r2 ' 

(3.46) 







.2/ 2g r 







G'L 



tx 



f , 5\ G[ 



L2 + [2g 2/ + r 



tx 



+ 



L2 ' e2 r2 ^ 



r2^ 



Note that Ax = and = C, a constant, automatically solve the equations of motion. 
At the linearized level, this solution is equivalent to a diffeomorphism x — > x + Ct. From 



the definition (3.29) we see that this diffeomorphism also acts on the boundary metric. In 



order to retain a Minkowski metric at the boundary, we will set the boundary value of Gtx 



to zero. Another solution to eq. (3.46) is 

Ax = ^C , G 



tx 



-f{r)L^/r^ + 1)C 



(3.47) 



with G a constant. This solution corresponds to the coordinate transformation x ^ x — Gt, 
t ^ t — Gx, an infinitesimal boost along the x direction, under which the boundary theory 
Minkowski metric is invariant. Note that this solution satisfies Ax^vq) = 0. There are two 



more linearly independent solutions to eq. (3.46). One of them diverges at the horizon. 



located at r = rg, and we obtain the remaining one numerically in section |4} 

With the boundary conditions we have discussed, the near boundary expansions of Gtx 
and Ax take the form 



G, 



tx 



+ 0{r- 



Ax, 



je 



2^3 



+ 0{r- 



(3.48) 



2^4 ' V ^ ' ^ 2r2 

Recall that we are working in a gauge in which the scalar is real. If we act on the series 



expansion (3.48) with a gauge transformation of the form 

Ax -^Ax + dx i^x) , V ^ ^e^"^" 



(3.49) 



we will be in a frame where the boundary value of the gauge field vanishes. In this frame we 



can use (3.44) to identify ^ with the gradient of the phase of the scalar. We will see shortly 



that j is equal to the spatial component of the U{1) current and t gives us the space-time 
component of the boundary theory stress-energy tensor. 

As before, we can compute the partition function for this configuration. Using the 



on-shell action (3.25) with the bulk action (3.26) and boundary action (3.38), we find 



lnZ = (P-072) 



Volume 



(3.50) 
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where we have kept only terms which are quadratic in Gtx and Computing the one- 



point functions from eq. (3.30), we find that eq. (3.40) receives corrections 



(3.51) 



See appendixIBlfor details. Following the strategy presented at the beginning of this section, 



once we have a solution to eqs. (3.46), we can can extract t, j and ^ from them and then 
use eqs. (3.43a) and (3.43b) together with ^ = VxP to obtain pg. 



4 Numerical Results 



In order to compute first, second and fourth sound explicitly, we solved eqs. (3.32) and 



(3.46) numerically. To obtain the isotropic solution (3.32), we found it convenient to use a 



shooting algorithm. We specified g'^rg), ip^ro) and 0'(ro) and looked for a solution where 
either the 0{r^^/'^) term or the ©(r-s/S) 

term in a near boundary series expansion of 
Tp vanished. Generically, for such a solution one will find /(r) = /or^/L^ + 0{r) near the 



boundary. To obtain the near boundary behavior for / as specified by (3.33 ) we rescaled the 
time coordinate such that t — > tjyfjQ- This rescaling amounts to the shifts /(r) f{i^)/fo 
and 4){r) — > (t>{r)/y/jQ. To solve for the linear perturbations of and Gtx a shooting 
algorithm is not needed; changing A(ro) and Gtxiro) amounts to varying Vx and Ux but not 

Ps- 



To read off the various coefficients in (3.33) and (3.48) we compared our numerics with 
a near boundary series expansion of the bulk fields expanded to 0{r^^^). The comparison 
was carried out on a surface of constant r = 10" (near the boundary) with n an integer. 
For each solution we made sure that our fit changed by less than 0.01% under the rescaling 
n — > n + 1. To reduce numerical error we used at least thirty digits of working precision 
in our computations and n > 4. As the temperature was reduced both n and the working 
precision were increased in order to retain an accurate result. The main restriction on our 
numerical computations was CPU time. 



With the near boundary values of /, g, (f>, ip, Gtx and Ax in hand, we used eqs. (3.41), 



(3.42) and (3.48) to compute P, /i, p, ps, s and T. We could then compute second and 



fourth sound from eqs. (2.18) and (2.21). Some typical results are shown in figure [T| 

From figure |2] we find that for the A = 3/2 condensate, there is a critical value of q, 
qc{u), above which second sound is monotonic. For q < qdu), second sound exhibits a 
maximum. It appears that qc{u) increases with u. The behavior of fourth sound seems to 
follow a similar trend: For q > qc{u) it asymptotes, as expected, to first sound at small 
temperatures but for q < qc{u) it falls short of c\ at small T. A typical example of such 
behavior can be seen in figure [3j 

The A = 5/2 condensate also seems to have a critical charge qc{u) below which second 
sound has a distinct maximum. However, for relatively large q > qc{u) and small n, the 



An intermediate result in this computation is the second order contribution to Sbuik: 



j(2) 



1 

2Z 



- AxA^ 



r 

2^ 



./' 



Gtx , 3 GtxG't 
L4 ^ 2nl f 



+ 



1 



" Ax 



Gtx 
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A=5/2, q=9, u=0 A=5/2, q=6, u=5 A=3/2, q= 1 , u=3 




Figure 1: (Color online) First, second and fourth sound for a holographic superfluid given 
by the bulk action (3.26) and a scalar potential (3.27). From bottom to top, the blue line 
corresponds to second sound, the green line to fourth sound, and the red line corresponds 
to first sound computed directly from (2.8). Note that conformal invariance implies that 
first sound squared must be 1/3 for our 3+1 dimensional system. 



A=3/2, u=0 




q 




-19. 






-15. 


-12. 


0.30 


-9.0 


0.25 


-5.0 


0.20 


-3.0 


0.15 


-2.0 


0.10 


-1.0 




0.05 



A=3/2, u=l 



T/Tc 




q 




-20. 


C 




-15. 




-8.0 


0.30 


-6.0 


0.25 


-4.0 


0.20 


-3.0 


0.15 


-2.0 


0.10 


-1.0 




0.05 



A=3/2, u=3 



0.2 0.4 0.6 0.8 



T/Tc 




q 

-20. 
-8.0 
-6.0 
-5.0 
-3.0 
-2.0 
-1.0 



0.2 0.4 0.6 0.8 



T/Tc 




u 




-10. 




-6.0 


C2 


-5.0 


0.30 


-4.0 


0.25 


-3.0 


0.20 


-2.0 


0.15 


-1.0 


0.10 


- 




0.05 



T/Tc 



A=3/2, q=6 




-20. 
-10. 
-6.0 
-3.0 
-2.0 
-1.0 
- 



0.2 0.4 0.6 0.8 



T/Tc 



A=3/2, q=20 




-10. 
-5.0 
-3.0 
-1.0 



0.2 0.4 0.6 08 1.0 



T/Tc 



Figure 2: (Color online) Second sound for a condensate of conformal dimension 3/2, whose 



dynamics follow from the bulk action (3.26) and a scalar potential (3.27). 
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Figure 3: 
charge q ■ 



(Color online) Fourth sound for a condensate of conformal dimension 3/2 and 
= 1, whose dynamics follow from the bulk action (3.26) and a scalar potential 



second sound curves develop a local minimum. This is depicted in figure |4} It is tempting 
to speculate about the possible relevance of roton like excitations in connection with this 
minimum. The behavior of fourth sound for A = 5/2 condensates is similar to that of the 
A = 3/2 condensates. 

5 Discussion 

In figures |2] |3] and |4] we presented our main results for second and fourth sound for a holo- 



graphic superfluid with a bulk action as in (3.25). The prominent features of these results 



are the appearance of a critical value of the scalar charge qdu) and the low temperature 
Hmit of second sound which is at odds with Landau's prediction for an incompressible, non 
relativistic superfluid (see eq. ( |l.l[ ))p| 

The appearance of a critical value for the scalar charge may be related to the two 
mechanisms that can drive the phase transition to the superfluid state |10| [36l H2l H3] . In 
the first, more obvious mechanism, the coupling of the scalar to the gauge field via the 
covariant derivative of the scalar generates an effective mass term of the form 



m, 



cff 



+ q^g^^Al'ilP' . (5.52) 



Once the charge of the black hole increases beyond a certain critical value, the effective 
mass of the scalar will become low enough to generate an unstable mode and the scalar will 
condense. The second mechanism, noted in the current context in ref. |36j (but see also 
|42j ) ■ is associated with a violation of the Breitenlohner-Freedman bound |44| H5] . As the 
temperature is decreased, the black hole approaches extremality. The near horizon metric 
of an extremal black hole resembles AdS2 x 5^, and this AdS2 is only stable for scalars 
above a modified Breitenlohner-Freedman bound. Neutral scalars with a mass m satisfying 
—(f' /A < m?L'^ < —d{d + l)/4 will condense at low enough temperatures. It is possible that 



'^The existence of a critical charge was also suggested in ref. [37]. 
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Figure 4: (Color online) Second sound for a condensate of conformal dimension 5/2, whose 



dynamics follow from the bulk action (3.26) and a scalar potential (3.27). 



the critical charge we observe is the charge below which the second instability becomes most 
important in causing the phase transition. Another perhaps related possibility, investigated 
in ref. P6], is that qc{u) is associated with nonconformal IR behavior of the theory. 

In order to understand better the discrepancy between our low temperature results for 
c\ and Landau's prediction, we decompose into three pieces: 



3s 



2 2 



(5.53) 



where = T{ds / dT) ^ is the heat capacity at constant chemical potential and 



sT/{sT + /xpn) which, following eq. (2.19), we call the quasi-particle velocity. Note that 



eq. ( |5.53 ) is a trivial rewriting of eq. (2.22). As discussed in section 2.1 if the low lying 
quasi-particle excitations are phonons, one should find that = 3s and that ci 



9- 



Moreover, we expect from eq. (2.24) that the low temperature limit of fourth sound is first 



sound. Instead, as can be seen from figures |5] |6] and |3] in our model is not always equal 
to 3s, Cl cannot be identified with Cq, and C4 does not always approach ci. It would be 
interesting to see if degrees of freedom on the field theory side with a dispersion relation 
u = Cqk can be more directly identified. 

We note that when q is very large approaches 3s and Cg approaches unity for non 
vanishing temperatures. This behavior is expected; in the probe limit, where the metric 
does not backreact on the scalar field and gauge field, the spacetime is a neutral black hole 
in AdSs and the chemical potential and density are supressed by a factor of q^^ relative to 
the entropy and temperature. 

The behavior of fourth sound depicted in figure [3] can be understood in terms of the 
behavior of ps and p, cl. eq. (2.21). We find that for all values of q and u, p p,^ at 



low temperatures. However, as can be seen from figure [7] ps does not approach p for small 
q, at least not for the temperature ranges we've reached. Thus, following the arguments 
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Figure 5: (Color online) Plots of the ratio of the heat capacity at constant chemical potential 
Cfj_ to 3s for A = 5/2. The q = 1 curves for u = 0, u = 1 and u = 3 have not been displayed 
— they are of order 3 with a maximum around T/Tc ~ 0.1. The value of 3s/C^ for the 
A = 3/2 condensate is very similar to the one in the plot. 
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Figure 6: (Color online) Plots of = sT/{sT + npn)- According to [T^, for a relativistic 
phonon g 3jS, Cq IS the phonon velocity. 
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Figure 7: (Color online) The relative superfluid density for various condensates. Note the 
similarity to the fourth sound curves depicted in figure [3j 



in section 2.2 it is not surprising that fourth sound does not asymptote to 1/3. In these 



cases, second sound seems to vanish at low temperature. Similar behavior of second sound 
has been noted in mixtures of ^He and ^He where pn remains nonzero at low temperatures 
[281 [29]. 
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A An aside on two-point functions 

In section [2] we defined ps as the variable conjugate to ^. This definition is slightly different 
from the one in the superfluid literature (see for example |47J. In what follows we show 
that these two definitions coincide. 

Consider a liquid in which the superfluid component moves with a small velocity Vx 
while the normal component moves with a small velocity u^. In such a configuration, the 



energy momentum tensor and current are given by (3.43). Another way of defining pn and 
Ps is through a zero frequency and long wavelength limit of the Green's functions involving 
T^* and . Consider applying a spatially varying external velocity field of the form UiC^^^. 
Such a velocity field is canonically conjugate to the momentum density T**, i = 1, . . . , d and 
we model the effect by adding a term to the Hamiltonian 



5H= / d^xT^\t,x)Uie 



ikx 



(A.54) 



Time dependent perturbation theory determines the first order response of the system 



(A.55) 



17 



where the Fourier transform of the retarded Green's function is defined as 



(A.56) 



By rotational invariance, the Green's function decomposes into two components which give 
the response of the system to transverse and longitudinal waves 



Qti,tj 



(0,k) 



(k) . 



(A.57) 



When U is parallel to k, the driving force induces a longitudinal excitation and we expect 
both the normal and superfluid components to couple to it. In the case where U is perpen- 
dicular to k and sufficiently small, we expect that no momentum will be transferred to the 
superfluid component due to its vanishing viscosity. Thus, only the normal component will 
be dragged. The transverse wave can be generated by an experiment where an open ended 



cylinder is filled with superfluid, and is dragged parallel to its own axis. Using eq. (3.43a I 
we find that 



lim Gi^(k) 



lim G2'(k) = -{sT + iMp) 



k-+0 



(A.58) 



These expressions are similar to the corresponding non-relativistic expressions found in the 
superfluid literature. In the non-relativistic limit, the sT term is dropped because it is 



negligible compared to ^p^ and /ip. See eqs. (2.14) in the main text. 



A similar analysis can be performed to extract the same limiting behavior of the current- 
current two-point function. Instead of a velocity field, we couple the charge current to an 
external gauge field Ai = —nVie^^^. Instead of eq. (A. 54) we use 



6H = p d^xJ\t,x)Vie 



ikx 



(A.59) 



Performing a decomposition of the retarded current-current two point function into parallel 



and perpendicular components, similar to eq. (A.57), we find 



(A.60) 



If V is parallel to k, then the external field A is gauge equivalent to zero, and we expect 
no response from the system. On the other hand, if V is perpendicular to k, then V can 
be identified with the superfiuid velocity from before and we can read off the two-point 



function from eq. (3.43b), 



lim G^^(k) 
k^o 



lim^G;^"(k) = . 



k-»0 



(A.61) 



B Counterterms and one-point functions 



In this section we compute the one-point functions in (3.30) using the general prescription 
initiated in refs. [21 13] and elaborated on in refs. |38l 
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To vary the on-shell action we consider an arbitrary linear perturbation of the solution 



to the equations of motion. Under a shift by 6Aa, Sjab, and 6il^, the action (3.25) will shift 



by an infinitesimal amount. This amount can be written as a boundary term which we 
regulate by evaluating it at large r = r^o- We find 



(B.62a) 



for A = 5/2 and 



6S 



7 ( ■^UaF"'' 5Ab - l^n^'daS^ - 2,i>5^l:/L+ 



2Ki ' 



(B.62b) 



for A = 3/2. 



Using the background solution (3.33) and (3.48), we can express this result as a power 



series in 1/r. For example, examining the coefficient of 6 Ah, and working to linear order in 
Ax and Gtx, we find that 



(B.63) 



Inserting all such expressions into eqs. (B.62) and carrying out the functional differen- 



tiation in eqs. (3.30), we find, at the end of the day, that 

(r") = (3P + c,V'(°V('^) , 



'tX\ 



(B.64) 



(T^^) = {Tyy) = {T'') = {p- ciV^°V^^^) , 

(O3/2) = -2V(°) , (O5/2) = 2^(1) , 

with C3/2 = 5/4 and C5/2 = —3/4. Setting either ip^^') or il^^^^ to zero, we recover the results 
(3.40) and ( 3.51[ ) in the body of the paper. The trace of the stress tensor is 



Tf^ = _4ci^(0)V)(i) = (4 - A)(Oa)V'^*^ , (B.65) 
in agreement with the conformal Ward identity (see for example eq. (4.21) of ref. [40j ) . 
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